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Abstract. Motivated by recent experimental work by Folkenberg et al. we 
consider the effect of weak disorder in the air-hole lattice of small-core photonic 
crystal fibers. We find that the broken symmetry leads to higher-order modes 
which have generic intensity distributions resembling those found in standard 
fibers with elliptical cores. This explains why recently reported experimental 
higher-order mode profiles appear very different from those calculated numerically 
for ideal photonic crystal fibers with inversion and six-fold rotational symmetry. 
The splitting of the four higher-order modes into two groups fully correlates with 
the observation that these modes have different cut-offs. 
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The advances in fabrication techniques and structural control of both index- 
guiding as well as photonic band-gap photonic crystal fibers (PCFs) is believed to 
improve the agreement between experiments and theory based on ideal structures (for 
a recent review of the field we refer to Ref. [J and references therein). The recent 
investigation of micro-bending induced attenuation in PCFs |2] is an example of the 
successful merging of experiments and simulations as well as recent experiments jS] 
on multi-mode cut-off in small-core PCFs and the good agreement with theoretical 
predictions On the other hand, in the latter experiments the observed higher- 

order mode profiles were very different from those of an ideal structure despite the fact 
that the studied fibers appeared to have highly regular air-hole lattices. However, the 
relative strength of disorder may still be significant and in this Letter we numerically 
study the effect of such variations on both fields and propagation constants. 

We consider the type of PCF first studied by Knight et al. [HI for which the 
ideal structure has inversion and six-fold rotational symmetry, see Fig. ^ In the first 
row of Fig. 12 we show the electrical field intensity |£^^^p (the superscript emphasizes 
the ideal or "unperturbed" structure) of the first six (m = 1,2, 3... 6) eigenmodes 
at A = 780 nm of an ideal PCF with A = 1.4 /im and d/A = 0.67 corresponding 
to one of the fibers studied in Ref. 0. The results are based on a numerical fuUy- 
vectorial solution of Maxwell's equations in a plane- wave basis with periodic boundary 
conditions [7] and basis vectors which support the inversion and six-fold rotational 
symmetry of the dielectric structure. For the refractive index n we use n = 1 in air 
and for silica we use a Sellmeier expression with n = 1.45367. The fundamental modes 
(m = 1 and 2) have close-to-Gaussian profiles [HI E] • 

Historically, their degeneracy has been widely debated (see Refs. ^3 E] and 
references therein) but recently it was finally proved by Steel et al. jl(J| how their 
degeneracy follows from group-theory because E'f'' and e'^^ do not support the full 
symmetry. Recently, the degeneracies of higher-order modes were addressed by Guobin 
et al. P2|, but only in the case of a PCF supporting 10 eigenmodes. The general picture 




Figure 1. Schematic of the cross-section of a photonic crystal fiber. The dashed 
circles indicate small random displacements and variations in diameter of the six 
inner air holes. 
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Figure 2. Plot of the electrical field intensity. Row (a) shows results for the ideal 
structure and rows (b)-(e) are for different random disorder in the inner ring of 
the air-hole lattice. 



is quite complicated and the degeneracies and order depend on the number of guided 
eigenmodes |13| . In a picture of true eigenmodes this can be understood as follows; 
for a given number M of guided eigenmodes the eigenvalues are minimized with the 
constraint that the eigenfields are orthonormal. When the number is changed to, say, 
M + 1 the orthogonalization affects all M + 1 eigenvalues and eigenfields and thus 
the picture of degeneracies may in general also change compared to the picture for M 
eigenmodes. 

In our case we numerically find that while E'f' and E^^ support the full 
symmetry e''^^ and e''^^ do not. Group theory thus predicts that the PCF have 
a set of modes which are either pair- wise two- fold degenerate {m = 1,2 and m = 5,6) 
or non-degenerate (m = 3 and m = 4). This contrasts the situation of a step-index 
fiber (SIF) in the weakly-guiding approximation where the first higher-order modes 
are four-fold degenerate. The 3rd and 4th modes are nearly degenerate in the sense 
that our numerical studies show that the splitting pf^ — P^'' is much smaller than 
the spacings — and P^^ — P^^ to the two neighboring groups of degenerate 
modes. The grouping of the four higher-order modes into two groups is experimentally 
reflected in different cut-off wavelengths of the two groups ^ . Compared to the picture 



Mortensen et al, Small- core photonic crystal fibers 



4 



IDA V* 

# Y i" \ 



o 











-0.01 -0.005 



V' 



0\ 



o 



Y 

•V 
o 

P V 

Yo 





□ 


m=l 




m=2 


o 


m=3 


« 


m=4 


V 


m=5 


4 


m=6 







O 

0.005 0.01 



[27r/A 



Figure 3. Plot of distribution P(/3m) for the first six modes m = 1, 2, 3, 4, 5, and 
6. The dashed Unes show Gaussian fits. 



in the weakly-guiding approximation, the splitting of the higher-order modes in PCFs 
originates from the inversion and six-fold rotational symmetry and the magnitude of 
the splitting is driven by the high index contrast between air and silica. We note that 
for a standard fibre, a sufficiently high index contrast will eventually also split the 
higher-order modes. 

From perturbation theory (for the problem of shifting boundaries see Ref. \14\ ) 
it is obvious that for guided modes the six inner air-holes are of most importance 
to the optical properties whereas disorder in the air-hole lattice further out in the 
cladding mainly affects the guided modes in terms of e.g. leakage loss |15j . Thus, in 
the search of over-all qualitative effects of disorder it is sufficient to only let the six 
inner air-holes (see Fig. have small random displacements Ar^ (i — 1, 2, ... 6) with 
respect to the ideal positions of the air-holes as well as small random deviations Adi in 
their diameters. We consider the situation where the displacements and deviations in 
diameter are Gaussianly distributed (zero mean and widths cta and ad) and mutually 
uncorrelated, i.e. 

(An) = (Ad,) = (Ar.Ad,) = 

{Ar.Arj) = (7^% , (Ad.Adj) ^ ajS,^. (1) 
In the following we study a large ensemble (more than 1000 members) of fibers 
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Figure 4. Near field images at A = 780 nm recorded at the output of 2.0 meters 
of a PCF. Panels (a)-(c) shows the fundamental and higher-order modes. Panel 
(d) shows a microscope image of the cross-section of the PCF with A ~ 1.41 fim 
and d/X ~ 0.66. 



numerically. 

Since the structural imperfection is typically a small perturbation of the ideal 
structure we have (/3m) — f^m^ ■ The grouping discussed above thus also holds 
approximately for the average propagation constants (/3m) as we also have confirmed 
numerically (not shown). In Fig. |31 we show the distribution of the propagation 
constant for the first six eigenmodes at A = 780 nm of a PCF with A = 1.4 /im 
and d/A = 0.67. For the disorder we have used cta = = 0.02 x A corresponding to 
2% structural variations of the six inner air holes. The fundamental mode only has 
a small fraction of the field near to the air-hole silica interfaces compared to higher- 
order modes, see row (a) in Fig. |21 From perturbation theory one thus expects the 
fundamental mode to be less sensitive to disorder compared to higher-order modes. 
Indeed, from Fig. O it is seen that the distribution P(/3m) has a common width for 
m = 1,2 which is significantly smaller than the common width for the higher-order 
modes (m = 3, 4, 5, 6). 

For a non-zero a, symmetries are broken and symmetry-related degeneracies 
are in principle lifted. For m = 1,2 this splitting is often referred to as form- 
birefringence quantified by Ani2 — (/3i — /32)/fc where k is the free-space wave number. 

However, since (^Ef'^\5e\Ef^'^ ~ (£^2°'' |^£|-^2*''') fo'^ ^ scalar perturbation 5e there is 
no significant form-birefringence (Ani2 ~ 0) to lowest order in perturbation theory. 
Indeed, for the studied disorder strength we numerically find an extremely narrow 
distribution. The width is vanishing (within our numerical accuracy) corresponding 
to a delta-function distribution P(/3i — (32) — 5(/3i — (32). 

The effect of disorder on the fields is illustrated in Fig. |21 Rows (b)-(e) show 
examples of the electrical field intensity for four random configurations of the inner ring 
of air-holes. In general, we find that the modes in the ideal structure, row (a), deform 
into profiles having two maxima with a node in between similarly to the deformation of 
the modes in elliptical core SIFs ■ This picture correlates perfectly with the modes 
observed experimentally in Ref. [S]. Fig. 0] shows near-fields recorded at A = 780 nm 
for a PCF with A ~ 1.41 /xm and d/X ~ 0.66. Here, the fundamental mode, panel (a), 
corresponds to the linear combination E ^ Ei + E2 and for the higher-order modes, 
panels (b) and (c) corresponds to E ^ E3 + E4 and E ~ E5 + Eq , respectively. 

In conclusion we have demonstrated how even weak disorder in the air-hole lattice 
leads to deformation of the higher-order modes in small-core photonic crystal fibers. 
In the presence of disorder the higher-order modes resemble those in elliptical core 
standard fibers. Our findings explain the recently reported higher-order mode profiles 
and the grouping of higher-order modes correlates with the observation that they have 
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different cut-offs. 
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